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Motivation uncraka Pt 


¢ Unmitigated RFI (Radio Frequency Interference) 
can cause errors in science measurements 
— L- and C-Band: soil moisture measurements over land 


— L-, C- and X-band: ocean salinity, sea surface 
temperature, wind speed direction 


— K band: water vapor, liquid water 
¢ Approach 
— RF front end development for 18 GHz (K band) 


¢ These allocations are known to be corrupted by direct 
broadcast services 
— Digital back end to allow sophisticated RFI detection 
and mitigation techniques 
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SMAP TA H-pol filtered 


SMAP (Soil Moisture Active 
Passive) algorithms developed 
previously under ESTO (Earth 
Science Technology Office) 
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RFI from Geosynchronous Satellites +»... 
Reflecting from the Surface ci 
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Picture from David W. Draper, [1] 


The 18 GHz Channel sees significant RFI from surface reflections around 
CONUS (Continental United States) and Hawaii 
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Real Signal Kurtosis (RSK) ea 


Given a complex baseband signal z(n) = /(n) + jQ(n), the fourth 
standardized moment is computed independenily for both the real 
and imaginary vectors, | and Q, as was used in SMAP{[3]. 


[G-E[I])*] 


E 
Sa (E[(-E[I])])? 


_ _ E((Q-EIQ)4) _ 
5 BSKo = Gq-eiann? 


The test statistic, RSK [2,3] (Real Signal Kurtosis), is then defined as 


__ |RSKi|+|RSKg| 
2 


RSK 
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Complex signal kurtosis (CSK) [4,5] is used to improve ability of the digital 
radiometer to detect RFI. It makes use of additional information in complex signals. 


Given a complex baseband signal z(n) = I(n) + jQ(n), moments az, Of z(n) are defined as 
t= E|(z —E[z])*(z- E[z])*"| ,f€,m ERZ0 


With 0? = a,, , Standardized moments ep» can then be found as 
Aem 


22m ~ bem 
oO 
Leading to the CSK (Complex Signal Kurtosis) RFI test statistic used [4]. 
2 
02.2 —2- 22:0 


Cr = 
1 
1 + 5le20l 
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Independent Component Analysis vin 


e ICA [6] uses higher order statistics to perform blind source separation 


¢ This suggests it may be useful for separating RFI from Gaussian noise in the 
radiometry context, studied in [7]. 


e We assume noise and RFI are statistically independent sources, mixing is linear, 
sources are non Gaussian 


e Mixture model: x = As, observe x 
° S=Wx, Sis the estimated independent source 


Scieauuctors Observation vector x Estimated vector ‘s 


$1(1) 
S1(n) 
$,(n) 
S2(n) 
7 $3(n) 
Original S3(n) Estimated 
sources/signals S4(n) sources/signals 


S4(N) 
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ICA RFI Detection set 


ICA Output 


i max{ABS(RSK - 3)} 


Detector 
Output 
Step 1: Take Kurtosis of each estimated 
independent component vector Step 2: Select the kurtosis value that 

deviated the furthest from 3 
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ROC Curves and AUC ee 
¢ Each point on an ROC curve can be 


represented by the set {FAR, PD} ROC curves 


— {False alarm Rate , Probability of Detection} 


¢ ROC curves will generate from (0,0) to 


(1,1) by varying the threshold § 
oO 
¢ Poor detectors are close to the 1:1 line 3 
' ra) 
e Better detectors show higher PD and > 
smaller FAR 3 
rs Worse Detection 
¢ Figure of Merit = Area Under Curve (AUC) 2 : 
=-===Kurtosis Detection Algorithm - Fullband 
— 0.5<AUC <1 : * | <== Pulse Detection Algorithm 
— When AUC =0.5 detector does not work — Se eee — - = 
— When AUC = 1 the detector works perfectly . False Alarm Rate 


ROC curve example, from [8]. 


AUC=1 | gE | AUC=0.5 
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ICA Performance, CW d = 100%, N = 9000 More ICA results in [7], 
l l l 
{ —@- direct | RSK L generally a marginal 
ae improvement in detection is 
~~ fastica | RSK 
0.9 - @ -fastica | CSK L seen 
“O- robustica | RSK 
“fl robustica | CSK ICA Performance, CW d = 100%, N = 9000 
0.3 nccfastica | RSK 
—&* necfastica | CSK 
O ~O~ erbm | RSK 
= 0.7 ~O~ erbm | CSK = 
~O- cqamsym | RSK 
Xi, ~ Bl cgamsym | CSK 
0.6 S| ~O>~ cerbm | RSK = 
~O~ cerbm | CSK 
0.5 
0.4 
-4 -6 -8 -10 =12 -14 : ; ; 
INR INR at AUC = 0.75 
Various ICA algorithms are tested [9,10,11,12,13,14,15,16,17]. RSK = Real Signal Kurtosis 
No ICA pre-processing is done on ‘direct’ data sets. CSK = Complex Signal Kurtosis 
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Eigenvalue Approach ea 


¢ Two objectives: 


— Detection: Identify power measurements that 
have been contaminated with interference 
¢ The Minimum Maximum Eigenvalue (MME) 
approach, adapted from the cognitive radio context 
[10], is applied here for RFI detection in passive 
remote sensing. 
— Excision: Accurately guess what the power 
measurement would have been if the interfere 
were not there 
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Conceptual Signal Model oo vesu: 


Hypothesis Test / Signal Model 


eens Polar Orbit ee 
Stationary Orbit Radiometer 


: e 
sNR=—= P= E[(rIkI)?I 


w[k] ~N(0,02) =Thermal Noise 
r[k] = RFI 


(depiction not to scale) 
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IN MARYLAND 


Given our sampled signal x, The Eigenvalues of the 


covariance matrix are found 
x(n) =x[nM+i-1] i=1,2,...,M 


A,>A.>°>dh 
x(n) = [xy(n), x9(n), xy (n)]" _—- a 


SHE OGD. The test statistic is then formed 


Ry, = E[R(n)X" (n)] 


L-24+N, 


1 
Ry(Ne)=— ) R(n)K"(H) 
° n=L-1 


n 
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Estimate oe 


Scale the minimum eigenvalue of the covariance matrix to estimate 
the variance of the Gaussian thermal noise. The limiting distributions 


from [19] help derive the scaling factor. 


R.A, > Ay > nD Amn 
a Amin = 0°(1 — Jy) 
iim, Ane =O ay) 
_— N 


a2. = J... ——__>- 
W m" (/N, — VML) 
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Wideband RFI —5 QPSK 
Channels 
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MIME Detection Results earmennirs 


Compare detection ali ed 0000, M=4,L=4 


—— MME CW 
—©— MME 1CH QPSK 

i MME 5CH QPSK 

—— Kurtosis CW 

—©— Kurtosis 1CH QPSK 

~~ Kurtosis 5CH QPSK 

—&— Spectral Kurtosis CW 

—— Spectral Kurtosis 1CH QPSK 
—<+- Spectral Kurtosis 5CH QPSK 


Eigenvalue Detection method 
greatly outperforms all other 
methods tested (Kurtosis[2,3] 
and Spectral Kurtosis[20]) 


0.6 
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Sample Variance 
Compared to 
Eigenvalue Variance Estimate 


Eigenvalue Variance Estimate 

outperforms sample variance 

at interference levels of -14db 
INR and greater. 


Eigenvalue variance estimate 
accuracy depends on the 
complexity of the RFI 
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RFI Excision Performance ow: 


H1: MSE Comparisons, N=10000,M=4,L=4 


—+— MSE Sample CW 

—©— MSE Sample 1CH QPSK 

—%— MSE Sample 5CH QPSK Sample 
aie Variance 
—©— MSE Eig 1CH QPSK 

—3— MSE Eig 5CH QPSK 
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